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SUMMARY 


Formulas for the wave drag of shrouded symmetrical airfoils and 
shrouded bodies of revolution of arbitrary shape are derived by means of 
linearized theory. In the case of the airfoils the shroud consists of 
flat plates, and for the bodies of revolution the shroud is a cylindrical 
shell. The results obtained hold for a Mach number range dependent on 
the geometry of the configuration. Expressions are also given for deter- 
mining a class of body shapes for which the wave drag is theoretically 
zero. 


INTRODUCTION 


A body moving at supersonic speeds has a wave drag which can be cal- 
culated either from integrations based upon the pressure at the surface 
of the body or by means of a momentum balance over a control surface sur- 
rounding the body. The control- surface approach Bhows more clearly that 
the wave drag is related to the transport of momentum in the Mach waves 
created by the body. This approach also suggests the scheme of reducing 
or destroying the wave drag through the use of a shroud as first shown by 
Ferrari (ref. l) . With- such a shroud the waves are caught and reflected 
to the body surfaces where they may be absorbed without further reflec- 
tion. From the standpoint of the pressure exerted on the body itself, 
it follows that the reflected waves may strike the rear portion of the 
body in such a way as to provide a buoyancy to overcome the resistance 
of the body alone. The detrimental effect of the additional friction drag 
due to a shroud is not included in the present study. 

The principal object of the present investigation is to derive for- 
mulas for the wave drag of shrouded symmetrical airfoils and shrouded 
bodies of revolution of arbitrary sfiape. The airfoil is shrouded by flat 
plates and the body of revolution is shrouded by a cylindrical shell. 
Although many configurations are possible, the analysis here considers 
the particular arrangement where the shroud extends at least far enough 
forward to catch the Mach wave emanating from the body nose, and far 
enough rearward to cast Mach waves on the base of the body. As a special 
application of the results obtained, a class of body shapes, similar to 
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those given by Busemann (ref. 2) and Ferri (ref. 3 ) , are found for which 
the wave drag is theoretically zero. 

For either the case of the airfoil or the body of revolution, the 
analysis is based on the assumption of linearized theory and on opera- 
tional methods using Laplace transforms. Ward (refs, k and 5 ) has shown 
in considerable detail how operational calculus can be employed to treat 
bodies of revolution and quasi-cylindrical tubes. A similar approach will 
be followed here, except that Heaviside notation is not used. 

To the order of the analysis employed, discontinuities in the slope 
of the airfoil are admitted, while for the body of revolution the gradient 
of cross-sectional area and its derivative are assumed to undergo no 
abrupt changes. In actual practice, discontinuities producing fixed com- 
pression waves would certainly upset the accuracy of the results even more 
than for the airfoil or the body alone since the opposing surfaces offer 
the possibility of shock-wave and boundary- layer interaction. 


SXMBOLS 


»o 


A,B 


c d 


C D 


D 

f 

F 


B 

h 


speed of sound in free stream 

cross-sectional area of body in dimensionless terms. 



section drag coefficient. 



body drag coefficient, — — 

q 0 s f 

P - Po 

pressure coefficient-, — 


wave drag 

F 

function defined in equation ( 5 ), — — 

function defining upper surface of airfoil or generating 
curve of the body of revolution 

function satisfying relation (53) 

distance of shroud from x axis 
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j 1?Kq,Ki 


} 


Bessel functions (ref. 6) 


l 

m 


M 0 

P 


Po 


*0 

R 

S(x) 

S f 

to 

Uo 

x 

y 


Yi 


Po 


4 

I 

P 0 

<p(M) 


length of body, or chord length of airfoil 

23 0 * 

l 


Mach number in the free stream 
local pressure 
pressure in free stream 


free-stream dynamic pressure. 


P0 U 0 2 

2 


radius of body 
cross-sectional area of body 
frontal area of body 
maximum thickness of body 
free-stream velocity 

Cartesian coordinate in free-stream direction 

Cartesian coordinate, measuring vertical distance for airfoil 
and radial distance for body 

Bessel function (ref. 6) 

■s/Mo 2 - 1 

dimensionless variable introduced in equation (4), ^ 

Mach angle, arc sin r£- 


dimensionless variable introduced in equation (4), 


x 


free-stream density 

perturbation velocity potential in dimensionless terms, — 
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® ( x ,y) perturbation velocity potential 

ft influence function defined in equation ( 35 ) 

( ) Laplace transform of function 

MATHEMATICAL STATEMENT OF PROBLEM 



(Boh =1/4) 
Sketch (a) 



(l/4<Boh<l/2) 



Consider a symmetrical airfoil 
or slender body of revolution placed 
at zero angle of attack at Mach number 
Mo = U 0 / a o > 1, where a Q is the 
speed of sound in the free stream and 
where U 0 , the free- stream velocity, 
is alined with the x axis. The y 
axis measures vertical distance in 
the case of the airfoil and radial 
distance in the case of the body. 

The nose of the body is at the origin 
of the coordinate system. 

In the configurations to be con- 
sidered the airfoil is shrouded by 
two flat plates and the body of revo- 
lution is shrouded by a cylindrical 
tube. The shrouding plates (or tubes) 
are placed so that the distance h 
by which the shroud is removed from the 
x axis is such that l/k < p 0 h < l/ 2, 

where 0 O = cot p. = n/Mo 2 - 1 and l 
is the body length. The shroud is 
required to extend at least from 
x = (3oh to x = l - Bob, but it would 
produce no additional effects on the 
body if it were longer. Sketches (a), 
(b), and (c) show the geometry of the 
configuration with either the airfoil 
or body of revolution in three typi- 
cal arrangements. If Bob = all 

the Mach waves from the forward por- 
tion of the body are reflected onto 
the rearward, and when Bob > l/2 
there is no effect of the shroud since 
the waves are reflected behind the 
body. 


(Boh >1/2) 

Sketch (c) 
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The upper surface (or generating curve) of the airfoil or "body is 
assumed given by the function 

y = F(x) (1) 

It will also be assumed that the body closes at both ends (i.e., that 
F(0) = F(l) = 0) and that the thickness-length ratio of the body is suf- 
ficiently small relative to the Mach angle (j. that linearized theory 
applies. As a result, the perturbation velocity potential $(x,y) 
satisfies the partial differential equation 



2 

8 $ _ v 8$ _ 0 
dy 2 y dy 


together with the boundary conditions 


( 2 ) 



where the parameter v equals 0 and 1, respectively, for the airfoil 
and body of revolution. One then has the problem of finding a solution 
of $, and from this to determine the drag of the configuration. 

If dimensionless variables g, tj, cp defined by the relations 

g = x/p 0 h , tj = y/h , cp( g,n) = <&U,y)/u 0 h (!)■) 

are now introduced, the expression 
• for the body surface, as given by 
equation (l), becomes 



with m = 2f3 0 h/l, and sketch (b) 
becomes sketch (d) . The differen- 
tial equation is then 

(1/2 < m < I) 
Sketch (d) 
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86 2 ^ 2 


v 8cp _ 
T) dT) 


0 


with boundary conditions 


(6) 




(7) 


Operational methods based upon the Laplace transform are suited for 
treating the basic differential equation ( 6 ) for either the airfoil or 
the body of revolution. Denoting the Laplace transform of a function 
g(fi,n) by g( s ;il) Where 


00 

I(s;ti)=/ e" s ^g(|,Ti)d| (8) 

\Jo 

and employing this transformation for relations ( 6 ) and ( 7 ), one obtains 
the differential equation 


s 2 q5 


d 2 ® v d® 
dri 2 “ ^ ^ 


= 0 


together with the boundary conditions 


(9) 


and 



^^vPof 



f v df 



( 10 a) 


(10b) 


Once <p has been determined from equation ( 9 ) , the drag of the shrouded 
figures can be calculated. In order to carry out the calculations, it is 
convenient to treat the airfoil and body-of -revolution problems separately. 
The analysis for the airfoil offers little in the way of novelty but will 
be given first since it illustrates the framework of the methods employed. 
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AIRFOIL WITH SHROUDING PLATES 
Evaluation of Wave Drag 


The solution of equation (9) for the case of the airfoil (i.e., 
when v = 0) is 


<p(sjil) =a(s)e" ST l + b(s)e ST l (ll) 

Since vertical symmetry exists in the flow field, attention can he limi ted 
to the tipper half of the figure. Over the forward part of the airfoil the 
downstream inclination of the outgoing waves indicates, furthermore, that 
the second term in the answer (ll) can he deleted. Imposing the boundary 
condition (10a) , one then has 


t 

-saCsJe -81 !]^ = sf , a(s) = - f 

and it follows that the perturbation velocity potential satisfies the 
relations 


S(sjT]) = - fe~ ST l 

<p(M) = - f(| - tj ) 

or, in terms of the physical variables. 


(12a) 


<D(x,y) = - ^ F(x - p Q y) , 0 < x - p Q y < ml - 2)3^ (l2b) 

where m = 20 o h/l, l/2 < m < 1. 

The flow around the forward portion of the airfoil is given by equa- 
tions (12); in order to predict the flow around the rearward portion, 
however, it is necessary to determine the nature of the incoming waves 
from the plate. The velocity potential of these incoming waves can be 
obtained from the second term in the right member of equation (ll). 
Boundary condition (10b) requires that vertical velocity be zero at 
T) = 1. From equations (12) and (ll), therefore, one has 


sb(s)e s1 ^J T |_ 1 = - sfe" 87 !]^! ; b(s) = - ?e -28 
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Hence the potential cp 2 of the incoming waves is determined by 

$ 2 (s;ti) = - fe“ s ( 2_T ^ 

= - f(S + t \ - 2) 


(13) 


and the potential in the region of the plate is 


<p = - f(| - tj) - f(| + rj - 2) 


If the potential cp over the rear of the airfoil is written in the 

form 

= cp 2 (6,n)+ q>s(5>n) 

and qj 3 is assumed expressible by terms of the form a 3 (s)e -ST l, the 
boundary condition (10a) yields 





= sf 


Thus a 3 (s) = - fe -2s - f, so that 

<P 3 (M) = - f(S - T) - 2) - f(g - ri) 

The potential over the rear of the airfoil is then given as follows: 
q>(S,Tl) = - f(l ~ t\ - 2) - f(| + T] - 2) - f(| - Tj) (life) 

or finally 

$(x,y) = - J 2 [F(x - p Q y - ml) + F(x + Poy - ml ) +F(x- Poy)] , ml < x- 0^ < l 
Po 

(14b) 

It remains now to find the expression for wave drag on the surface. 
In the physical variables, the total wave drag is 
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w 


where g Q = p 0 U 0 2 / 2 is the dynamic pressure in the free stream of density 
p Q . With the substitution, from thin-airfoil theory, Cp = - (2/Uo) (d$/dx) , 
the pressure coefficient on the airfoil surface is 


°p 

C p = — [F*(x) + 2F»(x - m2)] 
Po 


0 < x < ml 1 


m2 < x < 2 
(1/2 < m < 1) 


J 

and the drag, expressed in coefficient form, can be written as 



( 16 ) 


If the reflecting plates of the configuration are designed for a fixed 
height h, formulas (l6) give the drag coefficient of a shrouded airfoil 

of arbitrary shape for Mach numbers Mo >->/l + 2 2 /l6h 2 . When 
0 o h = 2/4 (m = l/2), the waves are reflected from the forward portion 
of the airfoil onto the rearward as in sketch (a) , and for m > 1 the 
reflected waves do not intersect the airfoil. Formulas (l 6) thus yield 
results for an airfoil with or without shroud. 


Consider, as a simple example , the shrouded biconvex airfoil section 
whose upper surface is defined by 

2t Q 

F(x) = -yg- x(2-x), 0 < x < 2 (17) 

where t Q is the maximum thickness of the airfoil, 
given by formulas (l6), is then 


The drag coefficient. 
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„ _ 16 /I 
d " 3 Po V 

(W 3 - 6m + 3) , 

1/2 < m < i 


id ) 2 

m > 1 

c 

d 300 V 

l) ’ 


(18) 



with the parameter m = 2p Q h/i. It 
is seen in sketch (e) that for this 
special case, the value of for 

1/2 <1 < 1 is always less than its 
value for m > 1, and that a relative 
minimum occurs at m = l/'/lf. If the 
ratio - h /i is fixed, formulas (18) 
give the drag of the shrouded airfoil 
as a function of p Q for p Q > i/4h. 


Shrouded Airfoils Having Zero Wave Drag 


Determination of shape .- Although formula (16) for the drag coef- 
ficient will in general he greater than zero, there exist classes of 
airfoil shapes for which the drag is theoretically zero. Since the two 
integrals in the first expression in equations (16) can never he negative, 
the necessary and sufficient condition that the drag vanish is that 


F*(x) = - F* (x - mi) , ml < x < l 
F*(x) = 0 , (l - m)l < x < ml 


(19a) 


After integration, these relations become 


^ F(mi) - F(x - mi) = F(x) - i F(mi) , ml < x < l 


> (19b) 


F(x) = F(mi) , 


(l - m)i < x < mi 


Thus the airfoil can he drawn in an arbitrary manner from x = mi to 
x = i and the forward portion of the airfoil shape in the interval 
0 < x < (l - m)Z is determined; the portion in the interval 
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(l - m)Z < x < mZ is 
flat and equal to . 
F(mZ). Sketch (f) 
shows the geometrical 
construction of the 
profile. The upper 
half of such an air- 
foil is equivalent to 
the lower wing of the 
linearized version of 
a Busemann biplane 
arrangement, while the 
lower half is the upper 
wing. 



Sketch (f) 


For the special case when m = l/2, as in sketch (a), equations (l9b) 
can be written in the form 

F(x) - | F (i) = 2 F (l) " F ( x " l) * Z/2 < x< Z (19c) 

If, moreover, the airfoil is assumed to be symmetrical fore and aft, 

F(x) = F(Z - x) and equation (l9c) becomes 

F(x) - |F0 * |F0 - f (± - x) , 0 < x < Z/2 (I9d) 


In this event the forward half of the airfoil has odd symmetry about the 
ordinate of the quarter-chord position and, similarly, the rear half has 
odd symmetry about the three-quarter chord position. It is also found 
that the pressure distribution on the airfoil has fore-and-aft symmetry. 


Drag for off-design condition .- A shrouded airfoil whose geometry 
satisfies relation (19a), however, will only have zero drag for some 
particular value of the parameter m, say hiq. If such an airfoil is 
moving so that the parameter m is different from mo, formulas (l6) 
for the drag coefficient become 


c d 


c d 


c d 



l/2 < m < mo < l' 
l/2 < mo < m < 1 ) 


9 


m > lj 

(20) 
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As an example, let the upper 
surface of the airfoil (shown in 
sketch (g)) he defined by 


F ( x ) = " x ) > 

0 < x < 1/2 

= ^ (Z - x) 2 . 

1/2 <x <1 


J 


( 21 ) 


If the design of the configuration is such that 2poh = z/2 (i.e., 
hlq = l/2), relation (19a) is satisfied so that there is no drag. When. 
1/2 < 2p 0 h < l, however , the drag coefficient for the airfoil is given, 
from equation (20), by 


c d 



ifln 3 + 6m 2 - l) 


l/2 < m < i 

' 


( 22 ) 



m > 1 


A plot of PoCfl/Cto/l) 2 against m 
is shown in sketch (h) for values of 
m greater than l/2. Below m = l/2 
the calculations become more involved 
but the same general method is 
applicable . 


BODY OF REVOLUTION WITH SHROUDING TUBE 
Evaluation of Wave Drag 


The solution of equation ( 9 ) for the case of the body of revolution 
(i.e., when v = l) is 

q>(sji}) = a(s)Ko(sT]) +b(s)l 0 (sT]) ( 23 ) 

where Watson’s notation (ref. 6 ) for the Bessel functions Kq and I Q is 
used. Over the forward portion of the body the wave system is outgoing 
and the solution can be formulated from the first term in the right 
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member of equation (23). In order to impose the boundary condition (10a), 
one employs the relation 

= - sa(s)K!(sT)) - - 

OTJ ‘I 

the last expression holding by virtue of the fact that t) is swinil on 
the body surface. The boundary condition then yields, for f(0) = 0, 

a(s) = - 
- 2ir 

where the quantity A is 

a(i) - p 0 rf 2 (j) = = j«!(|) 

P Q h 2 P Q h 2 

S(x) being the cross-sectional area of the body, 
forward portion of the body is thus 

- - -gj s ^( 8 )K 0 (bt]) ( 2 6 ) 

From equation (2 6) and the boundary condition (lOb), the potential 
<p 2 of the incoming waves from the shrouding tube can be calculated. If 
the Laplace transform of q) 2 is assumed expressible in the form 

dm 

9 2 ( s ;ti) = b(s)l 0 (sTi) , - sb( s) I jl( stj) 

Oq 

the boundary condition at tj = 1 yields 


(24) 


= S(x). 
Po^ 


(25) 


The solution over the 


sb(s)l!(s) 


s 2 5 .(s)K 1 (s) 

2it 


3 


Thus 9 2 is given as follows: 


- b(s) = - 


sA(s)K 1 (s) 

2jtli(s ) 


9 2 (s;n) = - 


sA(s)K 1 (s)I 0 (sti) 

2jtli(s) 


&Pp _ _ s^A( s)K]_( s)l j^( st]) 

^ 2nli( s) 

(27) 


T 
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Let the potential cp over the rear of the body he written as 
cp(g,Ti) = <Pa(S^) + cPsCi^) where q> 3 is assumed to he given hy the 
terms of the form a 3 (s)Ko(sTi) . In the previous case of the airfoil, 
the normal gradient of cp 2 at the airfoil surface was of the same order 
as the Imposed boundary condition and it was necessary to take proper 
regard of tp 2 when 'equation (10a) was satisfied. In the present case 
the normal gradient of cp 2 is of higher order (at the body surface) in 
comparison with the contribution of cp , and the boundary condition is 

o 

satisfied within the accuracy of the theory by the relation 

a 3 (s) = - ± sA(s) , A(g) = *p 0 f 2 (£) (28) 

The solution over the rear of the body is therefore 

_ = _ sA(s)K 1 (s)Iq(sti) _ sA(s)Kq(sti) ^ ^ 

9 2jtli(s) 2 it 

In terms of the physical variables, the drag D of the configura- 
tion can be written as an integral of pressure over the surface of the 
body: 


_D_ 

*o 


■I'm. 


2 *F(x) g dx 


(30) 


Since for slender bodies of revolution the pressure-velocity relation 
becomes 

P - Po 2 1_ 

q Q U 0 dx U 0 2 Vdyy 

it follows that the drag integral is, in terms of the dimensionless 
variables. 



where the Laplace transforms of (p are given by equation (26) for 
0 < | < 2 and by equation ( 29 ) for 2 < £ < 2/m < b. Proceeding to the 
evaluation of the drag integral, we find from equation ( 26 ) that for 
small values of tj 



MCA TN 3718 


15 


|| = - i s 2 A(s)Ko(stj) 


g2 A(s) fin M 

2jt \ 2 





Sip _ _1_ 
St) 2jt 


s 2 A(s)K 1 (sti) 


sA(s) 

2jtn 


J 


(32a) 


where y Q is Euler’s constant. From equation (29), similarly. 




M 

Stj 


f“ A(s)Ko(stj) - |~ A(s) Jalfi Io( s n) 


5 f + ? o) - g «•) 


A( s)Ki( ST]) - A(s) Io’(sti) 

2j( Il(s) 




~ sA(s) 
2kj\ 


(32b) 


From these expressions, the inverse transform can ‘he written explicitly 
on the surface of the body in the following form: 


for 0 < g < 2 

| - s ¥ - e r s J^ - El)aEl 

Sep _ 1 dA _ df j 

Stj 2jtp 0 f(g) dg dg J 

for 2 < g < 2/m < ^ 

| ■ s A " tE)l11 ¥ • S ' Sl)aEl • I A " (E 

pi - 

Ij A”*(g 1 - 2)a(g - gjdgi - I A"(0)n(g - 2) 
fk = M 

St) dg 


(33a) 


2) O 


> 

( 33 ^) 
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where fl(|) is an influence function whose Laplace transform is 




eS K].(s) 

e _8 Xi(s) 


1 


(34) 


and where the usual restrictions A(0) = A* (0) = 0 have now been placed 
upon A(|). 

The evaluation of the inverse transform of fl(s) is given in the 
appendix. It is shown that over the range of 0 < £ < 2 the relation 


fl(g) = - | In + T(g) 


( 35 ) 



holds, where T(g) is a power series 
that is convergent over the range. 
Sketch (i) shows fl(g), and the 
accompanying table lists values of 
the function and its component parts 
for 0 < g < 2. 


£ 

- 1 1 nlli 
* 2 

T(!) 

flU) 

0.20 

0.033538 

0.122638 

0.156176 

.1*0 

.071029 

.21(6073 

.317102 

.60 

.113534 

.1(03432 

.516966 

.80 

.162601 

.568452 

.731053 

1.00 

.220636 

.755683 

.976319 

1.20 

.291665 

.970676 

1.262341 

1.1(0 

.383237 

1.220196 

1.603433 

1.60 

.512300 

1.512421 

2.024721 

1.80 

.732936 

1.85711(8 

2.590084 

1.95 

1.174207 

2.157066 

3.331273 

2.00 

00 

2.266000 

00 


Equations (33) now must be sub- 
stituted into equation (31). The 


resultant expression is the drag of 
the configuration. The linearity of the term 8cp/8g in the integrand. 


moreover, permits one to write drag in the form D = where 

Dbtj is the drag of the body as it exists alone in its own induced flow 
field and is the drag of the body attributable to the induced 


effects of the shrouding tube. 
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Up to this point in the analysis no explicit use has been made of 
the assumption that the body closes .and, in fact, straightforward evalua- 
tion of leads to the expression given by Ward (ref. If-) and Frank! 

and Karpovich (ref. 7) for open-ended bodies with finite slopes at the 
base. If, however, the body closes or has zero slope at the end, there 
results the simpler drag formula of von Karman (ref. 8) . 


2 /m 2 /m 

vr - - £/ A "(£i)in|s - ix|a tl 


o 

2/m 2/m 

-5/ I A '( 6 )A-( gl ) |£-Mi 


o o 

The interference drag is then 

2 /m 




dA 

« d| 


= - h 


2 /m 

2 p n ^$2 dA 


-I 


2 —la d| 

dg ag 6 


(36) 


where dcp 2 /dg is given by the last two terms in the right member of equa- 
tion (33b). Substitution into the equation for the interference drag 
yields 


2 /m £ 

5^=h 2 ^ A* (|) jjA m (g - 2) + A”(0)o(| - 2) + J' A»"C| 1 -2)fl(| -l^dgijdg 


(37) 


so that the formula for the total wave drag for a closed body of arbitrary 
shape becomes 

2/m | 

~ = h Z J A‘(g)|A"(£- 2 )+A"( 0 )n(g- 2 )+f A”» ( £3. - 2 )fl( £ - £i)d£ij d£ 


2 

2 /m 2 /m 

J A* (g)d£ J d£i , 1/2 < m < 1 


hf 

2 St 


> 


o 
2 /m 


2/m 

’ £ = i/ A,( 6 ) «[f^ dU ' 


m > 1 


( 38 a) 
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In terms of the physical variables, using equation ( 25 ), one may write 
the drag in the final form 

l x-m Z ' 

|^=J's*(x)|s n (x-mZ) + S M, (0)fl^- S m (x- ml - xjfi^^dxajdx + 

S'(.x)Sxf 2^-dx x , l/2<m<X > 




( 38 b) 


where m = 2p 0 h/Z . If the shrouding tube of the configuration is designed 
for a fixed radius h, formulas ( 38 ) give the drag fo r a shrouded body of 

revolution of arbitrary shape for Mach numbers Mo > Jl + Z 2 /l 6 h 2 . When 
m = l /2 the waves from the forward - portion of the body are reflected onto 
the rearward as In sketch (a) . For m > 1 there iB no effect of the tube 
nnfl the formula for the drag is the same as that of the body without 
shroud. 


As an example, consider the body of revolution whose generating curve 
is defined by 


R(x) = — ^ x(Z - x) , 0 < x < Z 

“L 

where* to is the maximum thickness (see sketch ( i) ) • Here 



S(x) = x 2 (Z -x) 2 

, v .. 

S’ (x) =5 ^! x (Z-x)(Z-2x) 
“h 

S"(x) = ^SL. (Z 2 - 6xZ + 6x 2 ) 
Z 4 

ItSirt 2 

S’”(x) = ^2_ (2x-l) 

Z 4 j 


(39a) 


(39b) 


with frontal area Sj = jrt 0 2 /4. For m > 1 the second formula in 
equations ( 38 b) yields. ;• 
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Integration of the first relation in 
formulas (38b) for l/2 < m < 1 is 
carried out numerically and the 
results are shown in sketch (k) . 

For a fixed value of the ratio h /l, 
this sketch gives the drag of the 


configuration as a function of 
Po for Po > l / kh ‘ 





5 fi 7 .8 .9 ID 


m 

Sketch (k) 


Shrouded Bodies Having Zero Wave Drag 


In general, the value found from the drag formula (38) for a shrouded 
body ’trill be greater than zero. It is possible to obtain from this for- 
mula, however, an integral equation whose solution will yield a class of 
body shapes for which the wave drag is zero. 1 

Derivation of integral equation .- In the particular application to 
be considered here, we take 0 o h = if 4. (m = l/2) and assume that A"(0) = 0. 
A further assumption to be made is that the body has fore-and-aft symmetry. 
This implies that A(g) = A(4 - g) and, since no discontinuities are 
allowed in the slope of the meridian section. A* (2) =0. Under these 
conditions the first relation in equation (38a) becomes 


= h 2 ^A«(g)|A"(S-2)+ J'A"*(i 1 -2)fl(i-i x )d| 1 + iy' dg x Jdg 

(*L) 

A sufficient condition for the drag to vanish is that the expression 
within the brackets be zero, resulting in the relation 


A n (S-2 )+J A« , (g 1 -2)ft(g-g 1 )d li+lj' J~i~ d&1 = 0 > 


2 < g < 4 
(42a) 


The solution of this integral equation for the function A(g) is required 
in order to find the area distribution S(x) = p 0 h 2 A(g) of the body 
whose drag is canceled by the shrouding tube. 

1 While this part of the analysis was being carried out, Graham, 
Beane, and Licher (ref. 9) published a paper treating essentially the 

same problem from a different point of view. 
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Solution of integral equation .- Equation (42a) can also "be -written 
E(i) = A"(i) + J - i jaii , o < i < 2 (42b) 


where 


E(5) * - y/ (~ h i - + 1 - 1 


t(i) = ~ «(£) 

<3-5 


(43) 


In dealing with this integral equation it is advantageous to apply the 
calculus of Laplace transform because the equation contains the convolu- 
tion of A” and ijr. 

Taking the Laplace transform of both sides of relation (42b), one 
obtains 


or 


E(s) = s 2 a(s)[ 1 + sft(s)] 


E(s) 


s 2 A(s) 


^Ki(s) 

jte“ s Ii(s) 


(44a) 

(44b) 


For a homogeneous equation of the type (42b) it is usually convenient to 
assume that the unknown can be expressed in the form 


A"(|) 



7(i 2 )g(S - £ 2 )ai 2 


(45) 


where the functions g and y are to be determined. Equation (44b) then 
yields 


E(s) = 7 (s)g(s) 


^Ki(s) 

jre _s Ii(s) 


(46) 


Now the quantity in the denominator of this relation can be canceled 
by simply assuming 


7( s) = ite“ s I 1 (s) 


(47a) 
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or 


r(i 2 ) = j -} s - t '0 < | 2 < 2 
^l 2 (2-£ 2 ) 


= 0 , | 2 > 2 
so that equation (46) reduces to 

E(s) = g(s)e s K!(s) 


or 



1+ Si 

^ ii(2+ ii) 


g(i- ii)aii 


(47b) 


(48a) 

(48b) 


From equations ( 45) , (47b) , (48b), and the first relation in equa- 
tions (43) , there results the following integral equation for deter- 
mining the function g: 


f 7 ^-^.,= g(8-fix)as x 

>/6i(2+ix) 


= - i [*\— i 
itJ Li -2- 


ii l-2+l 




1 - i 2 
n/ | 2 (2 - la) 


g(ii- i 2 ) d i 2 


(49) 


After integration and some manipulation, relation (49) reduces to the 
homogeneous equation of the first kind 



0 < | < 2 


with unknown 

Gi(ii) = g(Sx) + g(2-Si) 

and kernel 


k(i,ii) 


_I_ +I rJi.. _ tan' 1 / (2 ~ Sx)(6- 

^S-Si)(2+6-6x) ^ (2 + i-ii)ix 


(50) 


(51) 


(52) 


This expression obviously has at least the trivial solution 
or 


Gx(lx) = 0, 
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S(£i) “ - g(2-6i) , °<Si<2 (53) 

requiring simply that the function g in equation (4-5) he odd about 
ii » 1. For a function that satisfies this condition, it follows that 
a formal solution of the integral equation (42) can he written as 

A"(g) = r g(g - gjdii , 0 < | < 2 (54a) 

or 

J, 

A«(g) = J n/£i( 2 - U)g(&- £i)a£i (54h) 


In order to avoid obtaining bodies with negative areas, however, the 
function g satisfying relation ( 53 ) must be chosen so that 

A (r) =~f [(£i-l) *Vli(2- lx) + cos" x (l“ £i)]g(l- li)d|i 

( 5 ^) 

is greater than or equal to zero. These formulas determine the shape of 
the shrouded body for the range 0 < g < 2 , while the rear half of the 
body is obtained by reflection about g = 2. There is of course an 
infinitude of functions for which the relation ( 53 ) holds, so that equa- 
tions (54-) actually furnish an infinite class of body shapes whose drag 
is canceled by a shrouding cylindrical shell. It will be shown later 
that this function g can be related to the pressure coefficient on the 
shroud. 

The quantities defined by formulas (54-) are continuous if g is 
piecewise continuous. It is apparent from these equations that the 
restriction A"(0) = A* ( 0 ) = A(o) = 0 is met, and it can be seen by use 
of equation ( 53 ) that the requirement A* (2) = 0 is also satisfied. 
Moreover, it can be shown that if A* is required to vanish -at g = 2, 
and if g is restricted to be piecewise continuous, the solution 
G'i(gi) = 0 used for equation ( 50 ) is the only one admissible here. 


A simple function that satisfies the functional relation ( 53 ) is 
g(ii) = ^ 2 li > 0 < gi < 1 

~ ^ 2 (li“ 2 ), l<gi <2 
Inserting this expression into the formula (54c) yields, for 1 < g <2, 


(55) 
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A(6) = Tjg jjsdl 2 - 8g + 5)cos _1 (l - £)+ (g-l)(2g 2 - 1 ^+ 15 ) n/|( 2- |) 


1}8(2- g)cos“ 1 (2- I) -16(| 2 -1 j.| + 6) */(! - 1)(3 - £)J 

and, for 0 < £ < 1, 

A (g) = -jjg ^3(^ 2 “8| + 5)cos" 1 (l-|) + (|-l)(2i 2 -lH+15) J&7S] 


(56a) 


( 56 b) 


A plot of the coor din ates of the body 3, 
is shown in sketch (z) for the spe- 
cial case when the maxi im -mi thickness 
ratio is 0.1 and the Mach number is 

</2. 

A simpler example is given by- 
taking 

g(ti)-Ad-Si) , 0<Si<2 - 3 I 

(57) 




\ 


x/l 


\ 


Tib 


Sketch (Z) 


Substitution into formula (54c) leads, for 0 < g < 2, to 

A(l) = j^3(- ^ 2 +i6|-13)cos- 1 (1- i)+ (- 2£ 3 + l4£ 2 -35l + 39) ^£(2- g)J 

( 58 ) 


a result which is the same as that obtained by Graham, Beane, and Licher 
(ref. 9). If the quantity 7 \ in equation (58) is chosen so that the 
body considered in reference 9 has the same volume as that of the previous 
example, the former body has a maximum cross-sectional area 0.9k times as 
great. 

Determination of body shape when a portion of it is specified .- For 
a function satisfying relation (53), it is seen that formulas (54) yield 
a. class of bodies having zero drag. Instead of prescribing the function 
g and then calculating A, however, it is possible to find the entire 
shape of the body when only a portion of it is given. This is accomplished 
by starting with equation (42a) . 

Equation (^t2a) can be put in an improved form by translating the 
axes so that the origin is at the body midpoint. To this end, set 
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-2 


-I 


g = cr + 2 , li=(Ji + 2 

A(g) = a( 4- g) =A(2 + cr) =A(2- a) =B(a) = B(-ct) 


BW 


L 


Sketch (m) 


(59) 


To be determined 
_ /-Specified 


As shown in sketch (m) equation (42a) 
can then he interpreted as an inte- 
gral equation in which B(cr) is 
specified in the range -2 < a < -1 or 
1 < a < 2, and it is required to find 
B in the range 0 < cr < 1 by inverting 
the expression 


>* (< ’ ) • - v/ 


2c 


— 5 2 

cr 2 - cr 1 


(60) 


where jjl(ct) is known and is given by 


ix(o)-B"(a-c/)+J’ B"(2- «!)♦(» - J (&) 


with ljr(cr) = dfl/dcr. A final change of variables (tf2= T ) reduces equa- 
tion (60) to the form 


yUr) = _ 1_ r 1 (d/dTj.) [2 nTFI( d/dr^B^) ] 

W7 " 2* J t " r x 


(62) 


The singular integral equation (62) is a familiar one in aerodynamics. 
Since 



B(^) dT 1 = B , (l) 


its inversion is 



kirJT 


1 

yr^L J 



or 
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<*> 

Equation (63) can be integrated with, respect to cr; by inter- 
changing the order of integration in the last term and imposing the con- 
dition B’( 0 ) = 0 , we find 

i 2 B«(l) . 1 T 1 / t 2 n/1- oZ-oJl- r 2 2 

B*(ct)= — — sin x 0 + = / ti(T 2 )ln ===== r = cLt 2 (64) 

* * J TjWl-0 2 +0s/l-To 2 


Integrating once again with respect to cr finally yields the expression 


B(o) -Bd)*. 25 '^ jVl- g 2 + g sin - 1 ff-ij 

1 p 1 f \ t t 2 ^1 - cr 2 - 0 ^ 1 - Tj 
i / h(t 2 ) cr In -5 - ===== 

r 2 S/1-0 2 + cr Vi - Tj 


/l-ff 2 - -j 1 - T g ; 


- T 2 In ^t 2 

s/1- cr 2 + si - t 2 2 


for the value of B in the interval 0 < a < 1 when it is prescribed for 
1 < 0 < 2. 

As an example, specify B”(0) to be 

B’^ff) =~h*J 0(2 - 0) (66) 

in the interval 1 < 0 < 2 and find its value for the range 0 < 0 < 1. 

It can be shown that 

0 0 

J B"(2 - 0iH(0 - 0i)d0i = 7> J s/ 0i(2 - 0i)t(0- 0i)d0! 


= - s/ 0(2 - a)J 


so that, from equation (6l) , 


M f 2 tol 

* J 0 2 - 0i 2 
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Substitution of this quantity into the inversion formula ( 63 ) leads, 
after integration, to the result 



or finally 



0 < 0 < 1 


( 70 a) 


Integration of this then gives 



Sketch (n) shows a plot of the coordinates of the body when the m a ximum 
thickness ratio is 0.1 and Mq = *f2. 



Sketch (n) Sketch (o) 

If the function B(cr) is specified in the range -1 < 0 < 0 or 
0 < 0 < 1, equation (l»2a) may also he employed to find B(d) for 
-2 < 0 < -1 or 1 < 0 < 2. (See sketch (o) .) In this case we may write 
the relation (lj 2 a) in the form 
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ff 1 

B”(2-c r) - r B»'(2-ff)n(ff-ff 1 )afl 1 +^ f B ' ,(tTl) 2 dg]_ — 

J o a - 

1 2 

J B«»(2-g 1 )ft(g-g 1 )dg 1 = J' toi i 1 <ff <2 


where the last term cm the left, as well as the term on the right is 
unknown. Applying the inversion formula for the airfoil equation under 
the condition that B"( 2 ) = 0 , yields the equation 


G(cr) = B"(g) + J' H(g,g 2 )B"(g 2 )dg 2 


where the known functions G and H are given /by 
r r 2cr Ik - g 2 f 2 hiCrg) / t 2 2 - 1 

— J ir^-i dT2 


hi( t 2 ) « B"( 2 - t 2 ) + B”(gi)t(T g - 2+g3.)dgi + ^ f dg x 

2-4 2 J 0 T * " I 


H(cr ,<ra) = t( t£ - 2 + ^ 


Relation (72) is a Fredholm integral equation of the second k in d and can 
he treated by well-known methods (see, e.g., ref. 10), although an inver- 
sion formula for the equation cannot , in general, be written in closed 
form as was done (eq. (62)) for the previous case. 

Drag for off -design condition .- If the body geometry satisfies the 
integral equation ( 42 a), the drag will be zero only at the design con- 
dition p Q = l/ 4 h (m = l/ 2 ) . The value of the drag of such a body for 
1/2 < m < 1, however, can be calculated by means of equation (38a) and 
will lie between zero and that for the body alone (m > l) . Substitution 
from equation ( 42 a) into formula (38b) gives the equation 
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Consider, as an example, the shrouded "body whose geometry is defined 
by equation (58), that is. 


s ”« [<3 - > f)if (l-f)- co^C 1 - 4 f)] ' 


S"(x) = S"(Z - x) , 


l/2 < x < Z 



For m = l/2 the configuration has 
zero drag, since A" or S" is a 
solution of the integral equation (lj2a) . 
Integration of equation (7*0 is carried 
out for this particular case for 
1/2 < m and the results are shown in 
sketch (p) . Calculations for m < l/2 
become exceedingly involved and cannot 
be obtained from the analysis, although 
the same general method is applicable. 


Pressure Distribution 


Pressure coefficient of shrouded body of arbitrary shape .- It is of 
interest to determine the pressure distribution of a body of revolution 
shrouded by a cylindrical shell. From equations ( 33 ) > the pressure 
coefficient 

_ P ~Pp _ g 1 /WY 2 2 dcp /dg/N 

■P *3-0 U 0 U 0 2 \$yy 0o dg \8 ti/ 
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on the same "body in the range 0 < £ < 2 is given hy 


°p ( 7 ?a) 


and, for 2 < § < 2/m < 4 , hy 


°p* 111 A A»(o) fl ( 8 - £ ) 

p i 

^a”(i-2 )+~J A"*( 6 1 - 2 )n(i-i x )ai 1 (75b) 


In terms of the physical variables, using equation (25), we write these 
formulas in the final form 

/ \ a P x 

Cp = - ±L_L*i In -( 53. ) +i ~ / S m (x 1 )ln(x- x^dxx ; 0 < x < ml 

* ^ 2 \dx/ it ^ X J — — 

(76a) 

C P " - ^ 311 ¥-(s) 2+ i s/ s»(x 1 )aii(*-x x )ax 1+ s»(o)i>^- a)+ 

x-mZ 

S"(x-mZ) + J* S ,n (x~mZ - Xx^^Yy dx x j mZ < x < Z (76b) 


In the interval 0 < x < mZ the pressure coefficient is the same as that 
of the body without shroud, while for ml < x < Z it is affected by the 
last three terms in formula (76b) . One sees from equations (76a) and 
(76b) that at x = mZ the coefficient has a jump equal to S"(o) , which 
means that the pressure is not continuous unless the body has a cusped 
nose. 

For the particular example of the body whose generating curve is the 
parabolic arc 

R(x) =~r x(z -x) , 0 < x < Z (77) 

b 

one obtains for the forward portion of the body 
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°p - <^) a [- 1 + 16 f - 22 <1- 6 f + 6 f)] . 0<x<ml 

(78a) 


and for the rearward portion 
v2 


C P = 4 (t) [- 1 + l6 f -22 Y5 + 23t(l + 6m+6m 2 + 6 ||-6 | -12 y m)- 
(x/l)-m 

12* / (■ | _ 2m - 1 - 2z)fl^| z)dzj ; ml < x < l 


(78b) 



-x/1 


A sketch, of the body of revolution 
having max-imum thickness ratio equal 
to 0.1 and a plot of the corresponding 
pressure distribution at Mach number 

■J2 are shown in sketch (q) for 
m =* l/2, and also for m > 1 (when the 
"body is without shroud) . The drag 
of the shrouded "body is shown in 
sketch (k) . 

Pressure coefficient of zero- 
drag body .- When the body geometry 
is such that relation (42a) for zero 
drag holds, formula (75b) for the 
pressure coefficient on the rearward 
portion of the "body (i.e., in the 
interval 2 < g < 4) reduces to 



df 

,d£ 


so that for t 0 < x < l / 2 

°p = - (, x ) In ^ - (^j S"( Xl )ln(x -xjdx! (79a) 
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and for 2/2 < x < 2 

°P = - ^ ^§“ - (§) " \ B“(x 1 )la(x 1 - x)axi (79b) 

x 

Since the hody here is symmetrical fore and aft, that is, since • 

S”(x) * S' r (2 -x), it can he seen from equations ( 79 a) and ( 79 b) that 
the pressure distribution on the body also has fore-and-aft sy mme try. 

It remains to relate the pressure coefficient on the shroud to the 
function g in formulas ( 54 ) . From the first relation in equation (32b) 


If ° ~ ' Mifij ' ^iC^KoC^J+^sJloCsrj)] 


so that on the shroud (i.e., at tj = l) 


Po Portias) 


__ sA(s) 

0 o *li(s) 

Taking the transform of equation (54a) , one obtains 

s 2 A(s) = Jte“ s I 1 (s)g(s) 
or 

g(s) = 

Thus from equations (80) and (8l) 




( 80 ) 


( 81 ) 


( 82 ) 


so that the pressure coefficient on the shroud is related to g by 


g(i) = eo0p’(i+s) 


( 83 a) 
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It follows from this equation and relation (53) that 


Po^p — J g(li)d-5i s 1 < £ < 2 
o 

3 -| 

PoCp* f g(Ei)«ili , 2 < £ <3 


Cpa+^^CpCS- £) 


that is, that the pressure distribution on the shroud has fore-and-aft 
symmetry about the body midpoint £ = 2. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif, Apr. 6, 1956 
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APPENDIX A 

EVALUATION OF AN INVERSE TRANSFORM 


The Laplace transform of the influence function 1 0 was expressed 
in equation ( 34 ) by 



e s Ki(s) 

.5te" s I 1 (s) 


1 


(Al) 


In this section, representations 


of the inverse transform 


fl(|) = lim 

to — 00 


1 

2*1 


a+ito 


J ' e z £ft(z)dz 


a-ito 


(A2) 


are given in two forms. The first form, which involves a series of 
Bessel functions, enables one to find the nature of the singularity. 
The other form gives the singularity explicity and is more convenient 
for computational purposes, as in the numerical integration of equa- 
tions (38), (63), and (76). 


Evaluation of the line integral (A2) can 
he obtained in the usual way by transforming 
the path of integration into a closed contour 
and applying the calculus of residues. The 
integrand possesses a simple pole at the 
origin and an infinite number of poles (roots 
of I x ) along the imaginary axis. The inte- 
grand also has a branch point at z = 0 due 
to Ki, so a closed circuit is chosen as 
indicated in sketch (r) . Since for a number 
a > 0 



Sketch (r) 


|ft(Re ie )| < ~ , -jt < 6 < ir 


(A3) 


The integral taken over the ABD and GHI of the circle c x of radius R 
can be shown to go to zero as R — > 00 . According to Cauchy’s residue 
theorem, the evaluation of equation (A2) is then 


1 This influence function was derived and first investigated in some ’ 
unpublished work of Max. A. Heaslet of Ames Aeronautical Laboratory. 
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n(i) = ^r m d) -Q(i) 


(A^) 


where 

Q(l)= 2 ^j- f e z £ft(z)dz + JL f e z £ft(z)dz + A_ 11m J e geei0 n(ee i0 )d(ee ie ) 
DE FG 

(A5) 


and r m is the residue of e z £fl(z) at z = ±iAm, the quantity A m being 
the mth root of J a (A) =0. We now consider the two terms on the right 
of equation (Ah) separately. 


Since 

A [zl^z)] = zl 0 (z) 
dz 

a particular member r m under the s umm ation sign in equation (a 4) is 
given by 


1 £+2) Kl (i7^) e" iAm( ^K^-jAm)' 

L lo(i^m) lo( 


m JtiAi 


L ] 


Use of relations 


KCtiAm) = ± If Yi( ^ ; I 0 (±±\i) = J 0 (*m) 


then yields the result 


00 

I-I 

m=i 


YxCAm) 


cos(g+ 2) Am 


(A6) 

(A7) 


(A8) 


for the total contribution of the residues. 

The last term of equation (Ak ) , which represents the contribution 
resulting from integration over the path DEFG in sketch (r) , may be 
written from equation (A5) in the form 
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Q(t) 


2n s l 


11m 
e -> o 


' T e gte lg K 1 (te^ ff ) 1 

I * L lit*® 1 *) "j 


dt + 


f e £te- i]I e 2te- i \ i ( 
J t Ix(te- ; 


ijt ) 


- it 


it 

dt + i/ e6 ee 

it 


10 


e 2 ee i9 Kl ( ge 10 ) ^ 


IiCee 10 ) 


30 


If the relation 


K 1 (te ± i) _ K x (t) ± 
I 1 (te ±i ) Ii(t) 


(A9) 


(AlO) 


is employed, equation (A9) then gives 


r co 

QU) = ±- lim &t- 2 / • 

^ e o{_ J 6 


»-t( i+ 2 ) 


dt + 


if 


k (8+2)eeW K : 


i( ge16 ) J 

iCee 3 ^) J 


or 


Q(fi)-| [* + |-(6+2) 2 + ln 2 (6 + 2)] 


(All) 


(A12) 


Finally, using relation (A4) , (A8), and (A12) , we have the equation 

0(6) -i [ic£ y cos(| + 2)V-it--|+U+2) 2 -ln 2(£+2)] (A13) 


for the inverse transform of the influence function. 

The nature of the singularities in this expression may he found 
heuristically as follows. Except for the last term, which has a loga- 
rithmic infinity at £ = -2, the only possible singularities in equa- 
tion (A13) come from the series involving Bessel functions. Since the 
singularities due to this series are not affected by omitting the first 
N terms, we substitute the asymptotic expansions (ref< 6) 
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a — .. — aL 1 1 — - 

8rt(m +^J 1283^(01+^ 


J + • 


sin(7v m -^)+ . . .J 

YlW * [“ COS ( Am "C ) + 8fc S K Am " ?) + * * *] 


> (A14) 


so that 

Yi(Am) 


■I 

m=N 




eMUt2fc -^ CO8(i C 2)>m [- 1+0 fe). 

m=N 


^ cos ^m+^(|+ 2 )jt ^ sin(m+^ 

w m + i “ m + i 

m=N I 4 . m=N Ij. 




when N is sufficiently large. Now 
00 sin 


+ i 1 + din 

y . ^ 4 

4 m+ i 2 


^ 2 1 - sin 

4 


(A15) 


(A16) 


which, near | = ±2 behaves like ln(2 + i)/(2- |) . One can see from equa- 
tion (A13) that the function 0 has only a logarithmic singularity at 
5=2.” 


EVALUATION OF ft(£) BI ALTERNATE METHOD 


The expression (A13) for the influence function ft is unsuitable 
for numerical work, so that it is desirable to have an equation which is 
more convenient and which also isolates the singular part. We therefore 
seek an evaluation in the form 


ft(|) = - | In ^i + T(|) , 0 < £ < 2 (A17) 
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where T is a convergent power series 

CO 

T(|)=^a j3 | j (ALB) 

i=o 

If the function 5(z) in the line integral 

a+iw 

fi(S)=~- J e z £ft(z)dz (A 19 ) 

a-iw 


is expanded as an asymptotic series } the desired series expansion of 
ft(g) can he obtained through term-by-term integration. 

For large values of |z|, asymptotic expansions of K 1} are 



IlU)s 7fe[ p (-0- le '^@] ' 

where 


J arg zj < 3*/2 
-it/2 < arg z < jt/2 


> (A20) 




m=x 


(2z) m 




1 *. (i,m) , 
l6 m (2m- l)(mi) 3 IT ^ x 


> 


J 


(A21) 


so that for small g equation (A19) yields 


a+ioo 
a- ioo 

The term e -2z F(l/z) in the denominator can be neglected if Re(z) 
(i.e., if a) is chosen sufficiently large and positive. Thus 


F 


g(l/z) 


F(-l/z) - ie -2z F(l/z) 


-1 


dz 


(A22) 
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a+ico 




a-ioo 


(A23) 


Expanding the "bracketed term as a series in l/2z, one obtains 


a+ioo 


n(s) “sr I eZl Z^ 


dz 

z m+i 


a-ioo m?=i 


or, after term-by-term integration 


tl(s) sr(£) > °- <5 < 2 

m=i 


(A2^) 


where the coefficients c m may be determined in succession from the 
relations 


ci=2bi 
c 2- Cxbx 

C3 = c 2 bx - Cxb 2 + 2b 3 


> (A25) 


= Om-xT 3 ! “ c m-2* 3 a + • • • + ('"l) Cxbm-i + [1 + ( _ l) + ]^m J 


Cm 


vm+ii 


If use is now made of the expansion 


ln(2- |) =]n 2 


00 -m 


(A2 6 ) 


m=i 


it follows that the series T in equations (A17) and (Al8) is given by 


T(t) = I = [o^jrtl 

m=i 


-» / \ m 
JD > 


0 < £ < 2 


(A27) 
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and that the influence function ft can he expressed in the final form 


a( S )=-iinl|i + 


r i [_2e lYif 

L m L(m - 1) I *JW- 9 


0 < i < 2 (A28) 
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